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Related Applications 



This application claims the benefit of Provisional Applications Serial Nos. 60/222,800, filed 
08/02/00, 60/242,232, filed 10/17/2000 and 60/288,273 filed 05/02/2001. This application is 
related to pending formal applications Serial No. 09/299,022 and 09/433,964. 

Previous filings by the author are included by reference for the entirety of their disclosures. The 
first is "Tilt-Compensated Interferometers," filed 4/26/1999, serial number 09/299,022. More are 
provisional applications serial number 60/107,060, filed 1 1/04/98, titled "FT-IR Signal 
Processing: Part I," serial number 60/1 19,429, filed 02/09/99, titled "FT-IR Signal Processing: 
Part II," a formal application entitled "Signal Processing for Interferometric Spectrometry" serial 
number 09/433,964 filed 11/04/99. Further provisional applications which are included for the 
entirety of their disclosures are titled "Interferometers and Interferometry," serial number 
60/228800, filed 08/02/00, titled "Interferometers and Interferometry: Part 2," serial number 
60/242232, filed 10/17/2000, and titled "Interferometers and Interferometry: Part 3," serial 
number 60/288,273 filed 05/02/2001 . The book by Griffiths and deHaseth, "Fourier transform 
spectrometry," ISBN 0-471-09902-3, is also included for the entirety of its content. 

Portions of the inventions disclosed here have been made under contract with the United States 
Federal Government through the Department of Defense under contracts DAAD13-01-P-0012 
and DAAD-05-96-P-3 102. The Government has certain rights in these inventions. 



It is an object of the present inventions to provide new interferometers, which are better than 
prior art in respect to stability, scan speed and cost of manufacture. It is an object of the present 
inventions to provide signal-processing techniqiies for interferometry, particularly very rapid 
scan interferometry. It is an object of the present inventions to improve the state-of-the-art in 
photometric accuracy of interferometric measurements. 



Background and Summary of the Invention 



Brief Description of the Drawings 



Figure 1 shows a diagram of a tilt-compensated interferometric spectrometer. 




Figure 2 shows a diagram of the reflector/beamsplitter assembly of Figure 1 . 

Figure 3 shows a diagram of a Michelson interferometer incorporating refractive path difference 
scanners. 

Figure 4 shows a diagram of a modified Michelson interferometer design incorporating refractive 
path difference scanners. 

Figure 5 shows a diagram of a bearing and support mechanism for producing a nutating motion 
of a refractive scanning element. 

Figure 6 shows a diagram of a remapping function applied to modify the apparent wavenumber 
spacing of the x-axis to actually wavenumber spacing. 

Figure 7 shows a diagram of a signal processing chain that has the effect of recording the inverse 
transfer function of a detector signal chain. 

Figure 8 shows a diagram of signals which are produced during a real-time measurement of 
detector saturation and correction for saturation. 

Figure 9 shows a diagram of an arrangement which may be used to measure and optionally 
stabilize the wavelength of a tunable solid-state laser to render it an effective wavelength 
standard for interferometry. 

Figure 10 shows a diagram of signal produced during measurement of the wavelength of a 
tunable solid state laser. 

Figure 1 1 shows a diagram of a permanently-aligned, tilt-compensated four-port interferometer. 

Figure 12 shows a diagram of a permanently-aligned, tilt-compensated four-port interferometer 
with both inputs derived from a common source, and the outputs recombined to effect optical 
subtraction. 



Figure 13 shows an alternate geometry to that of Figure 12 in which the functions of the separate 
beamsplitters and combiners are done by a single assembly. 

Figure 14 shows a diagram of a four-port interferometer according to the approach of Figure 13. 

Figure 1 5 shows a diagram of an optical path difference scanner that uses gas pressure to effect a 
variable optical delay. 

Figure 16 shows a diagram of an optical path difference scanner that uses an airbearing assembly 
to avoid shear alignment errors. 

Detailed Description 

Beamsplitter-roof reflector assembly 

Larsson has described in US Patent 5,650,848 a method for tilt-compensating the scanning 
mechanism of an interferometer. This prior art extends ideas described by Steel in the early 
1960's (see for example, "Interferometry," Cambridge University Press, 1967 W. H. Steel). A 
new approach substantially increases the throughput that is possible, while making it 
straightforward to produce very compact reflector assemblies and interferometers. The crux of 
the new invention is an assembly that consists of a beamsplitter rigidly mounted to a roof 
reflector. By roof reflector, what is meant is a pair of mirror facets which are mutually abutted 
and mutually perpendicular to interferometric tolerances. 

A slight disadvantage of the Michelson interferometer is that the two end mirrors in the arms are 
susceptible to misalignment with each other and with the beamsplitter. The Jamin (1856), 
Solomon (US Patent 5,675,412), Turner and Mould (US Patent 5,808,739), Woodruff (US Patent 
4,391,525) and several related designs provide an interferometer framework in which slight 
misalignments of the optical components are completely compensated. Thus, by adding one or 
more refractive scanning elements to such designs, it is possible to construct an interferometer 
that preserves intrinsic tilt-compensation and allows variation of path difference with minimum 
reduction of throughput. If the refractive scanning elements have faces that are parallel and flat, 



then they will usually not compromise the intrinsic tilt-compensation of an interferometer. The 
exception is that for interferometers that do not utilize path reversal, the shift caused by such 
optical scanning elements can cause shear, but not angular, misalignment of the beams in the two 
arms. 

Figure 1 shows a diagram of an interferometer which is permanently aligned and for which 
scanning of the optical path does not require interferometrically accurate rectilinear motion. 
Radiation from a source impinges on a beamsplitter coating which is supported on one side of a 
substrate. The beam reflected from the coating traverses to a hollow cube comer or other 
retroreflector where it is directed back through a portion of the beamsplitter substrate which is 
transparent. This traversal of the beamsplitter substrate provides compensation for the dispersion 
of the substrate material. The beam then propagates to a flat reflector which returns it on its path 
back to the beamsplitter. The radiation which passes through the beamsplitting coating traverses 
to a block which is rigidly affixed to the beamsplitter substrate. Preferably the block is a hollow 
roof reflector (two facets at 90 degrees to each other, both of these facets being oriented at 90 
degrees relative to the beamsplitting coating). If these two faces are arranged to be mutually 
perpendicular to the beamsplitting coating of the substrate, then three facets will form a cube 
comer, or lateral transfer retroreflector. The relevant surfaces of the block and can be coated with 
reflective materials by a variety of methods. The block is understood to be of interferometrically 
precise constmction such that after a reflection from the block the beam propagating to the flat 
reflector will have a wavevector that is exactly parallel to the beam from the retroreflector. Thus, 
the flat mirror will return the second beam back to the beamsplitter. It may be discerned that after 
recombination, both beams will have traversed an exactly equal thickness of the beamsplitter 
substrate. 

Two novel methods for fabricating the roof reflector and end mirror of the interferometer are 
described. These are wire electron-discharge machining (EDM) and fired ceramic. It is thought 
that by the choice of metals which have very uniform coefficients of expansion, together with 
wire EDM, that excellent precision substrates can be prepared. Preferred metals are aluminum, 
invar, and stainless steel EDM uses an electric arc to erode metal from a part. The electric arc 
produces very little damage to the part and particularly does not produce stresses in the parts that 
can cause adverse thermal distortion. Wire EDM can produce surface finishes that are excellent 
for replication. Replication is a known method for producing high-quality optical surfaces. 




Another preferred approach to preparation of substrates for the roof reflector is slipcasting of 
ceramic mixtures and molding of ceramic pastes (i.e., clays). By judicious choice of components, 
the thermal expansion properties of the ceramic can be tailored to be uniform in all directions 
and very small. In particular, the addition of ceramic microspheres is preferred because the 
resulting masses have uniform properties. Some ceramic compositions have particles that have 
non-uniform properties. If the particles are shaped such that they align preferentially (e.g., flakes 
or rods), then the macroscopic properties of the resulting masses may not be uniform. After 
casting or forming, it is understood that the resulting ceramic pieces are fired, then used as 
substrates for replication. 

This inexpensive interferometer is suitable for use with the nutating disk refractive scanning 
elements described below as well as a variety of other scanning methods. The beamsplitter with 
anti-parallel reflector assembly may be used as a component of the optical subtraction 
interferometers described below. The refractive scanning elements described below be inserted at 
any point in the optical paths in either arm of the interferometer. Motion of the cube comer 
provides a method for fine-tuning the zero path difference point of the interferometer, or it may 
be used as the optical path difference scanning element. Alternatively, rotating disks described in 
US Patent 5,898,495 may be interposed in the path between the beamsplitter and the cube comer 
reflector. It may be necessary to manufacture the final flat reflector in a stair step configuration 
to compensate the differences in path length of the two arms of the interferometer, thereby 
insuring that zero path difference is accessible by the scanning element or elements. Two views 
of the beamsplitter with roof reflector are also shown in Figure 2. These clarify the layout of the 
components on the beamsplitter. 

Nutating disk, refractive scanning interferometers 

The use of a nutating disk of transparent material as a refractive scanning element of an 
interferometer is disclosed herein. The principle of optical path difference scanning is that the 
disk is mounted such that a beam passes through it. During rotation, nutation occurs so that the 
angle of the disk with respect to the beam varies; consequently the optical path is varied. Such a 
disk must generally be supported by the edges so as to not block radiation from being transmitted 
through the disk. The disk is rotated approximately about the axis near where the center of the 
beam passes through it. The axis of rotation is not exactly perpendicular to the plane passing 




through the edge of the disk. Thus, the disk nutates or precesses as it rotates thereby causing the 
optical path through the disk to vary. The disk should have faces that are very flat and parallel 
such that after passing through the crystal the angle of propagation of the beam is unchanged. In 
cases where the beam is inverted and passed through the disk a second time, it may be acceptable 
for the faces to not be parallel, but they should still be flat to interferometric tolerances. 

A series of methods for mounting and spinning these disks are disclosed. Two preferred 
materials for construction of the disks are sodium chloride and potassium bromide, hereafter 
NaCl and KBr. The disadvantage of these materials is that they are relatively soft. However, 
calculations show that the deformation of KBr even when spun at 500 revolutions per second, 
suitable for measuring 1000 interferograms per second, will be quite small. Other materials are 
also contemplated for use as nutating refractive scanning elements. In particular, germanium is 
good because of its stiffness and high refractive index. Because of the high refractive index of 
germanium it should be used with antireflection coatings. Many other materials including silica 
(quartz), alumina, silicon, zinc sulfide and zinc selenide may be used depending on the desired 
spectral range, spin rate, retardation and cost range. 

The advantage of the nutating arrangement relative to other arrangements described in the 
literature (e.g., US Patent 5,898,495) is that the throughput can be much higher. In the case of 
Dybwad's design (US Patent 4,654,530) the substantial advantage of the nutating disk approach 
is improved duty cycle efficiency; however, the throughput is generally not significantly greater 
than that obtained by Dybwad. The throughput is significantly greater than that obtained by 
Davis (US Patent 5,883,713). The disadvantage of Davis' design is that the beam has to pass 
twice through a salt crystal, which is generally cubic. The duty cycle efficiency is consequently 
poor, as is the throughput. A point should be made that the rotation axis of the scanning element 
in Dybwad and Davis are both perpendicular to the scan axis described herein. In both the 
present case and US patent 5,898,495, the fiindamental scanning mechanism is nutation. Nutation 
of a transparent disk allows much higher throughput than nutation of a disk mirror. The reason 
for the difference is that with the disk mirror the beam must be reflected back to the disk mirror 
for a second reflection. The length of the return path is on the order of 3 d, where d is the beam 
diameter. In a nutating disk design, the total length in either arm of an interferometer must be on 
the order of 6 to 10 d. By using a nutating transparent disk, the total length in each arm of an 
interferometer can be reduced to the order of 3 d. Throughput can thus be approximately 10 




times higher because it scales as the square of the number of beam diameters of path traversed in 
each arm. 

Collins, et. al have disclosed in US Patent 3,955,083 the use of a nutating disk as an 
interferometric shaft angle encoder. The device differs considerably from the present disclosure, 
both in purpose and implementation. Collins described a disk mounted by a shaft through the 
center. Nutation of an optical element for the purpose of interferometric modulation has also 
been published by Crane in US Patent 4,743,1 14. Crane's art differs considerably from the 
present disclosure for several reasons. In particular, the motion of the beam produced by nutation 
of the optical element was not cancelled by a second pass in the reverse direction, but rather 
exploited as a means of scanning the field of view^. Further, Crane's element was a Fabry-Perot 
etalon rather than a disk for refractive scanning of optical path difference. 

A transparent nutating disk element can be used with a variety of basic interferometer designs, 
including the standard Michelson. Figure 3 shows the preferred embodiment with a Michelson 
interferometer in which one disk is located in each arm. The disks are oriented at approximately 
45 degrees relative to the beams. The disks are held by their edges in strong support cylinders, 
which may lend some precompression to enhance the stability of the crystal in a centrifugal field. 
The cylinders may be mounted in ball bearings for rapid rotation. The disk is oriented at a small 
angle relative to the rotation such that it appears to wobble or nutate during rotation. This has the 
effect of altering the path length for radiation passing through the crystal. A preferred 
embodiment has the cylinders supported by magnetic levitation. The same means can also drive 
very rapid rotation, by altering the magnetic fields at the support coils in synchrony. Either one 
or both disk elements can be rotated. If both are rotated, they may be rotated at the same rate, or 
at different rates. If one is rotated at twice the rate of the other, then interleaved scans of high and 
low resolution will be obtained. 

Figure 4 shows a diagram of a folded Michelson interferometer using a variation of an 
arrangement published by Steel many years ago. The advantage of passing the beam twice 
through the nutating disk is that any variation of angle caused by a wedge angle between the two 
faces of the nutating disk (i.e., non-parallelism) will be cancelled by the double pass. Thus, the 
nutating disk may be of less exacting manufacture. The disadvantage is that the extra length of 



the interferometer arms generally decreases throughput. Another disadvantage is that the disk 
must be twice as large to accommodate two beam diameters on its surface. 

If the angle between the nutating disk and the beam is approximately 45 degrees then the 
variation of path difference with angle is nearly linear for small nutation or rotation angles. It is 
found that for a 6 mm thick KBr window having a refractive index of approximately 1 .5, that the 
variation of path difference with angle for one disk with the beam passing in both directions is 
about 0.15 mm per degree of nutation (because the effect of 1 degree of nutation is -1 to +1 
degree of tilt). If a disk is mounted in each arm, then the maximum optical path difference (and 
resolution) may be doubled. 

If one disk is rotated at three times the rate of the other and is made appropriately thinner, then it 
may be used to linearize the change of optical path difference with respect to time. This idea may 
be furthered by adding two more disks and rotating them at 5 and 7 times the base rate to further 
linearize the change of optical path difference with respect to time. These disks would be 
correspondingly thinner. It is thought that the benefits of linearizing the scan are modest, given 
the signal processing approaches described herein and elsewhere. 

It is understood that there will be dispersion in the scanning element or elements. If a refractive 
scanning element is used in only one arm of the interferometer, then the interference signal may 
be strongly chirped. Mattson has published a method (US Patent 5,491,55 1) for signal processing 
which accounts for a strongly chirped signal. 

Figure 5 shows two methods for supporting and spinning an assembly that includes a nutating 
disk mirror. The figure shows the housing surrounded by a ball bearing. It is understood that the 
rotating assembly should be balanced by drilling material from the housing so that during rapid 
rotation vibrational forces are minimized. The disk should be housed in a suitable material. 
Carbon fiber composites are preferred because of their lightweight, high stiffness and small 
coefficient of thermal expansion. Many other materials may be used for the same purpose. 
Titanium is thought to be a good choice because of its excellent mechanical properties and 
lightweight. 



Signal processing for Dispersion 




A pending application by the author describes methods of signal processing which are also 
appropriate for use with refractively scanned interferometers. The general approach favored by 
the author is to record both the laser and infrared channels separately and simultaneously using 
suitable analog-to-digital converters or other technology such as that proposed by Mertz circa 
1965. These ADC's will be clocked from a common source so that they produce a set of matched 
simultaneous samples. An inverse transfer function may be applied to both channels to provide 
delay matching. The laser signal is then processed, preferably with a Hilbert transform or filter, 
to extract an accurate estimate of phase for each data point. The mirror position is proportional to 
this phase value. This set of data can be used together with the corresponding samples of the 
infrared channel to do either an unevenly-spaced Fourier transform (i.e., Lomb periodogram), or 
a spline or polynomial fit with subsequent conventional Fourier transform. Brault's method may 
also be used to produce excellent results. If the interferometer does not exhibit dispersion, the 
resulting spectral data are equivalent to those produced by conventional signal processing means. 
The signal processing steps detailed below will remove the effects of dispersion. 

In the case of a refractively scanned interferometer, spectral distortion will result from 
dispersion. All refractive elements exhibit some dispersion. The spectral distortion can be readily 
predicted and compensated given that the temperature of the refractive scanning elements and the 
relationship between dispersion and temperature are known. 

Doyle has found that it is not necessary to perform a correction for variation of temperature when 
using KBr as a refractive scanning element. Nevertheless, there are materials such as germanium 
for which it would be necessary. It is thought that for thermal and stability, the interferometer 
should be enclosed in a rugged metal housing. If the outside of the housing is even lightly 
insulated, the high thermal conductivity of the metal will tend to make the inside isothermal. In 
the preferred embodiment, the spectrometer has a small thermocouple or resistance temperature 
detector built into it so that the temperature of the metal shell can be continuously monitored. 
The thermal mass of the metal housing will insure that whatever temperature variations do occur 
are relatively slow so that corrections at each point in time will be accurate. 

The form of the error, in the spectrum, caused by dispersion is distortion of the x-axis, i.e., the 
wavenumber axis. If there were no dispersion, this axis would be as accurate as the reference 




laser stability. In order to adjust the spacing of the x-axis, it is necessary to fit a polynomial or 
spline to the data. The polynomial is then evaluated at unevenly spaced points along the x-axis 
which correspond to a grid of evenly spaced points in the absence of dispersion. This correction 
may be repeated for each temperature. The dispersion curve can be determined from literature 
coefficients, or from a NIST-traceable w^avenumber absorption standard. The traceable standard 
is preferred. 

Figure 6 diagrams the distortion of the axis and shows schematically the appropriate signal 
processing steps to remove the distortion. Variations of the peaks on the spectrum are grossly 
exaggerated for clarity. The remapping function can be readily calculated from the dispersion 
curve for the refractive scanning material at the operating temperature. The laser data give a very 
accurate indication of the thickness of the material. The modulation frequency in the top trace of 
Figure iS is on a linear scale and is referenced. As the refractive index of the material increases, 
the modulation frequency increases in exact proportion. To produce the bottom trace on Figure 6 
it is only necessary to modify the x-axis of the top trace by applying the remapping function. The 
remapping function simply adjusts the value of each frequency point by dividing by the refractive 
index at that point to produce an accurate wavenumber value. This leaves the x-axis data points 
unevenly spaced. If an evenly spaced wavenumber axis is desired, as is the convention in 
chemometric analysis, the next step is to interpolate to find an evenly spaced set of values. 

Signal Processing to Record and Invert Transfer Functions 

To effect the type of signal processing described above most effectively, the transfer functions of 
the sample and reference (generally, infrared and laser) channels must be recorded and inverted. 
A previous disclosure describes the use of LED's to probe the transfer functions. 

Figure 7 shows a powerful method of recording and simultaneously inverting the transfer 
functions of the sample and reference channels is to use two adaptive filters. The adaptive filters 
described here may be of the finite-impulse response type describe by Widrow and Steams 
(ISBN 0-13-004029-0). Each weight may be updated for each iteration of the filter computation 
cycle. The adaptive filter outputs are compared to delayed copies of the input signals driving the 
LEDs (assuming that the LEDs are operated in a region of linear response to their input signals). 
The error signal for each iteration may be computed from the difference between the delayed 




version of the test signal and the output of the filter. The error is applied as feedback to update 
the filter weight vector by a small amount at each iteration. Thus, each adaptive filter will 
converge to the inverse transfer function of the input channel, including the constant delay which 
is the same for both channels. The spectrometer should be controlled during the convergence 
time so that no spurious signals are introduced into the adaptive filter. Once the filter is fully 
converged, it thought to be sufficient to provide brief periods for updating the filter weight 
vectors to account for slow variation over time. These brief periods may occur between scans. 

Figure 7 details the approach. An LED or other emitter excites the detector. In a typical 
spectrometer, there would be both a reference laser channel and an infrared signal channel. Thus, 
the signal processing detailed in G would typically appear twice in an actual system. A waveform 
is applied to drive an LED. The output of the LED shines on the detector. The waveform may be 
calculated, as from a pseudorandom algorithm, or may be generated from a lookup table. A copy 
of the waveform that is applied to the LED is delayed by an amount sufficient that the 
corresponding output from the detector has reached the midpoint of the filter delay vector. The 
output of the filter is compared to the waveform that is applied to the LED (or a linearized 
version of the waveform applied to the LED). The reason that a linearized version may be needed 
is that the LED may have some non-linear response that can be precompensated by adjusting the 
input voltage/current to produce the desired output intensity waveform. In such as case, the input 
waveform differs slightly from the output of the LED. By comparing the output of the fiher to 
the delayed version of the input to the detector, the adaptive feedback which corrects the filter 
weight vector will force the filter to converge to compensate for the detector transfer function. 
The input waveform from the LED must not have frequency components that lie outside the 
response range of the detector. If this condition is not met, the filter convergence problem will be 
ill-posed and the resuhs unstable and inaccurate. This method would generally be applied 
between scans of an interferometer. During the time of scanning and data acquisition the filter 
weight vector would remain constant. Because the detector and preamp transfer functions would 
usually change only slowly over time, it would not be necessary to rapidly update the inverse 
transfer function unless adverse conditions are encountered. It is understood that the probing of 
the detector transfer functions would best be accomplished during times that the interferometer 
mirror is stopped between scans. 




For probing the degree of saturation of the detector continuously during scanning, a different 
approach is required. Detector saturation can be readily monitored by exciting the LED of Figure 
8 with a signal moderately higher in frequency than the highest frequency component expected in 
the signal or reference channel. In this case, the LED is exciting the detector simultaneously with 
the interferometric signal of interest. It is not thought that the reference channel generally 
requires compensation for saturation because the dynamic range of a sinusoidally modulated 
laser signal is quite small. In the composite signal, the higher frequency signal can then be 
demodulated to produce an estimate of the degree of saturation of the detector at each point in 
each scan. The magnitude of the recorded signal can then be modified for improved photometric 
accuracy by using the saturation information. To account for differences in the effect of 
saturation at different wavelengths multiple LED's can be used, each modulated with a different 
frequency so that the degree of saturation for a series of wavelengths can be simultaneously 
measured. 

Figure 8 shows a diagram of 4 signals obtained when exciting the signal detector with both the 
interferometric signal and the probe signal. The top signal designated A shows an interference 
signal which has suffered from the effects of detector saturation. The dotted line indicates the 
intensity that should have been recorded. The trace B indicates the interferogram with a high 
frequency probe modulation present. The trace C indicates the demodulated intensity of the high 
frequency probe modulation that clearly dips at the location of the highest intensity in the 
interferogram. The degree of the dip very accurately captures the degree of saturation. Trace D 
represents the intensity of the interferogram after division by trace C. This is the corrected 
interferogram which has improved photometric accuracy. 

LED's may be also be used to probe the transfer functions of the elements of an array detector. 
Silicon CCD arrays, mercury-cadmium telluride arrays, indium antimonide arrays, and a variety 
of other array detectors are used with spectrometers for imaging and other applications. 
Typically, these detectors produce an output proportional to the integrated intensity falling on 
them during a time interval. Thus, the time-dependence need not be probed. However, the non- 
ideal intensity response characteristics of these components may be probed and compensated on 
an element-by-element basis. These non-ideal characteristics include non-linearity and offset. 
The output of the LED's or other emitters may be homogenized by passage through a diffuser, 
integrating sphere, light pipe, or other optical element to insure that each element of the array is 




equally illuminated. The probe radiation for each element may be modulated independently as 
discussed below, but it is thought that use of equal illumination of the elements is a very practical 
method for probing non-uniformity. 

An LED or other emission source is controlled by a computer processor which may be a DSP 
processor. The intensity can be accurately controlled in stepwise increments. By adjusting the 
intensity through a wide range, and recording the output of each detector element for each 
intensity, the gain, offset and non-linearity of each pixel can be separately calculated. When 
spectroscopic data are recorded, the transfer function of each element may be removed from the 
recorded data by correcting for gain, offset and linearity of response using the recorded transfer 
functions of the detector pixels. The saturation of the detector can be probed by taking 
interleaved frames from the detector with the LED on in one, and off in the next. This is 
analogous to the use of phase modulation of the interferometer in step-scan imaging. 

The use of phase modulation in an infrared spectrometer may be used to reject multiplicative 
noise of the type associated with source intensity variation. This method requires that an image 
frame be recorded while the interferometer mirror is at one extreme of its motion. This signal is 
subtracted from the frame recorded with the mirror at the other extreme of motion. This would 
usually be repeated several times at each interferometer step and the results averaged. This 
practice will greatly reduce the effect of multiplicative signal variations of all types. 

LED's may be also be used to probe detector transfer functions of photoacoustic (PA) or 
photothermal (PT) detectors. PA and PT detectors are intended to be very sensitive to 
temperature variations caused by absorption of light. Unfortunately, they are also very sensitive 
to environmental disturbances as well as convection and other non-ideal characteristics of their 
components. LED's or other solid state emitters, either singly or in arrays, may be used to probe 
the instantaneous transfer fimction of the detector and each point within the detector sensitive 
area. Thus, changes in responsivity caused by temperature variation, aging, density fluctuations 
(convection) or other causes may be compensated in real time. The principal variation needed 
relative to the above discussion is that to probe different points in a photothermal or 
photoacoustic detection system simultaneously requires a multiplexed probe signal. For example, 
in the case of a photoacoustic sample that is mounted within a photoacoustic cell, a suitable 



approach according to the present invention is to modulate a probe beam which is impinged on 
the full area of the detector active area. 

The probe beam can be conveniently modulated by the use of a mirror array of the general type 
discussed by Messerschmidt (US Patent 6,034,370) and Tague (US Patent 5,923,036). 
Alternatively, a suitable probe beam can be generated by an array of LED elements that are 
assembled with hybrid circuit or integrated circuit technology. In either case, the probe beam is 
be imaged onto the active area. Each point in the probe beam is then modulated with its own 
characteristic modulation pattern that may be a sinusoid of a particular frequency, or may be a 
modulation described by a Hadamard matrix. By demodulating the photoacoustic detector signal 
for each of the patterns contained in the probe beam, it would be possible to discern the local 
amplitude and phase response of the sample-gas-microphone coupling system. This is thought to 
be useful, because, particularly in step-scan Fourier-transform photoacoustic measurements, the 
variation of response over time can be source of spectral noise. By compensating these variations 
during data acquisition, much higher quality spectra can be measured in a much shorter time. 

A source is collimated, then directed onto an array mirror device. The beam reflected from the 
array mirror device is directed into the photoacoustic sample chamber such that it illuminates the 
sample. By modulating each mirror pixel with a characteristic pattern, such as a sinusoid of 
known frequency and phase, it is possible to measure the amplitude and phase responses of the 
cell for each pixel independently. In the case of convection currents and hot spots that vary the 
temperature of the sample and the gas near it, it is thought that this information can be used to 
account for the variations in response. The composite signal from the detector, during a condition 
in which the response probing is active, but the input beam from a spectrometer is not active, is 
decomposed by Fourier transformation. As long as the probe signals are separated in frequency 
by more than the resolution of the Fourier transform, each probe amplitude and phase will be 
distinctly available. Ideally, all pixels are probed simultaneously so that information about the 
sample response is very dense. It is thought that the probing should be as quick as possible so 
that changes over short time scales can be compensated. 



Etalon Standard for Reference Lasers 



It would be very desirable to replace helium-neon lasers, which are typically used as frequency 
references for interferometry, with small, efficient semiconductor lasers. The power savings are 
up to 20 watts per laser replaced. This is a very important issue for portable interferometers. The 
waste heat generated by gas lasers has been determined to adversely affect the performance of 
some spectrometers. Unfortunately, for the most part, semiconductor lasers have a frequency 
stability on the order of parts per thousand, while the HeNe lasers have a frequency stability on 
the order of parts per million. The use of a small and inexpensive etalon (typically a quartz plate) 
to determine the absolute frequency of a laser is disclosed. An etalon with a 3cm"* mode spacing 
may be useful, but smaller or larger spacings can also be used. Using an estimated refractive 
index of 1.5 would indicate an etalon thickness of about 1 1 1 microns. The etalon may be used 
with or without an array detector. The etalon may be placed anywhere in proximity to the 
interferometer, in the path from the laser to the laser detector. 

Figure 9 shows a schematic of the optical and electronic layout appropriate for using such an 
etalon for fine control of the laser frequency. The source laser is controlled by a digital-to-analog 
converter which allows both a DC offset current to be applied, as well as an AC modulation 
current: The modulated output of the laser passes through an etalon and then through the 
interferometer. The optical output of the interferometer which has been modulated by the static 
position of the interferometer mirrors, as well as by passage through the etalon is then sampled 
by an analog-to-digital converter circuit. The resulting digitized signal is then demodulated by 
the DSP core which controls the overall spectrometer. The DSP calculates any corrections to the 
laser DC current to maintain a constant wavelength for the reference laser. The DSP also records 
the laser wavelength to use for correcting recorded interferometric data. 

Choosing the mode spacing of the etalon is a compromise between two conflicting criteria - the 
first being that the mode spacing be narrow enough that there would be significant intensity 
modulation when the wavelength is varied by current modulation, and the other being that the 
modes are sufficiently separated that there is no ambiguity as to which mode the laser begins in 
at startup. It is possible to use an independent reference to ascertain in which mode of the etalon 
the laser is operating. For example, the rotational lines of water are known to great precision and 
are often present in mid-infrared spectra. By calculating the position of the water lines, the laser 
wavelength can be very accurately determined. 




For semiconductor lasers, one powerful approach is to modulate the current in the laser to cause 
the wavelength to shift, and then decode the resulting intensity modulation to determine precisely 
the wavelength. A first order correction must first be applied to the data recorded from the output 
of the laser photodiode or photodiodes. The reason for the correction is that the intensity of the 
laser will vary approximately linearly with applied current. The intensity of the laser may be 
monitored with a second photodiode so that the laser intensity can be completely divided out of 
the values determined from the photodiodes. 

Many modulation approaches are possible - the modulation could be slow relative to the 
scanning modulation and done continuously during operation, or it could be done only between 
scans with a fairly high frequency. As the semiconductor laser current is modulated, the shifting 
wavelength will lead to an intensity modulation as the laser moves over one or more etalon 
modes. Different types of lasers will have different tuning ranges. If the current tuning can be 
done only over a small wavelength range (« 1 mode) then it is still possible to extract from the 
ratio of the fundamental and the first harmonic the exact location of the laser wavelength on the 
etalon modes. If the current tuning range is relatively large (on the order of 1 mode) then it is 
possible to drive it to an even or odd symmetry point on one of the modes and maintain it there 
continuously. 

The waveforms that are produced by wavelength modulation are Bessel functions. The frequency 
components of the various Bessel functions are very sensitive indications of the exact point in the 
etalon mode that the laser frequency is positioned. The amplitude of the Bessel functions is also 
a sensitive gauge of the wavelength range of the tuning. By comparing the amplitudes and phases 
of the various harmonics it is possible to determine the laser wavelength with excellent precision. 
This approach is capable of holding or determining the laser position to at least 1/1000 of a 
fringe, if not 1/3000 of a fringe or better. This approach gives a solid-state laser comparable 
frequency stability performance to a HeNe. The only cost incurred for each additional unit is the 
etalon. In the spectrometer design preferred by the author, the semiconductor laser is driven by a 
digital-to-analog channel from the DSP core, so no additional components would be required in 
the design. In the preferred design, the same DSP core is connected to the photodiode which 
detects the laser intensity as modulated by the interferometer. The additional frequency 
components which are produced by the various modulation schemes aimed at determining the 
laser frequency can be readily demodulated by the same DSP core. Since this DSP core, in the 




preferred design, also handle the interferometric data from the signal channel; any necessary 
corrections to the data can be readily recorded. 

Figure 10 shows the waveforms that may be expected at different points in the etalon modes. For 
example, there are four points of pure symmetry within one cycle of the mode spacing. These 
points are the positive and negative extrema (A and C, respectively), which have purely even 
symmetry. That is, the etalon modes viewed from this point are symmetric. The two points of 
purely odd symmetry are the zero crossings, B and D. 

Two important approaches which are thought to be complementary to the above description may 
also be used with lasers which cannot be conveniently wavelength modulated by current 
modulation. These are piezoelectric modulation of the etalon thickness, which will produce a 
similar effect. Another effective method to provide modulation for detection of wavelength is to 
use an electro-optic modulator at the output of the laser. It is thought that a standard HeNe which 
drifts by approximately a few ppm can be further improved for much less expense than the cost 
of a stabilized laser by the use of either of these methods. These methods are also applicable to 
stabilization or determination of absolute frequency of solid-state lasers of various types, 
including semiconductor lasers. With gas lasers, the initial frequency will be known to much 
higher precision than the semiconductor laser. Thus, the mode spacings of an etalon may be 
much narrower so that the frequency determination is correspondingly more accurate. 

The sensitivity to misalignment (i.e., non-perpendicularity of the etalon) is quite acceptable up to 
0.1 degree where it reaches about 15 ppm/degree. At 1 degree the sensitivity is 10 times larger 
and completely unacceptable. Thus, the etalon must be precision mounted with respect to the 
semiconductor laser. If the etalon is mounted in such a way that the angle can be modulated, then 
it is quite straightforward to modulating the etalon tilt to find a minimum in the thickness. The 
variation of refractive index and thickness with temperature can be readily corrected if the 
system temperature is probed as described herein. 

Optical nulling interferometry 

Optical subtraction, or nulling, of interferometric signals has been known for many years as 
practiced and published by Griffiths, Mattson (US Patent 4,999,010) and others. The main 
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advantage of optical subtraction is that the dynamic range of the measured interferogram is 
reduced. Dynamic range is particularly important in the near-infrared spectral region where the 
extra brightness of blackbody sources and the higher D* of detectors can make the dynamic 
range of interferograms considerably larger than 1 million:!. The performance of analog-to- 
digital converters was an important issue for mid-infrared spectral measurements for many years. 
Recently the performance of ADC's, particularly those of the delta sigma type, has reached the 
point that it can generally be disregarded for mid-IR measurements. However, the measurements 
made in the near-IR spectral region using interferometers can routinely achieve optical SNR 
levels well above 1 million: 1 up to the range of 1 billion: 1. Presently no ADCs can achieve this 
level of performance. It is impossible using current technology to digitize some interferometric 
NIR optical signals without seriously degrading the signal quality. 

The use of permanently-aligned, tilt-compensated four-port (two- input, two-output) 
interferometers to perform optical subtraction of signals is disclosed. The two inputs are also out 
of phase, so that if a suitable comparison source (or sample) can be found, then the difference 
between the inputs can be nulled. An example of a suitable standard would be a pharmaceutical 
tablet which is known to be good. Such a tablet would be a suitable standard for inspecting other 
pharmaceutical tablets for quality control reasons. By using a near infrared probe beam it would 
be possible to obtain the chemical signature of each tablet, thereby verifying the composition as 
it passed by on a conveyor-belt or was loaded into containers. Such quality control would be 
particularly useful in cases where more than one type of pill was produced at the production 
facility. For optical subtraction, the ideal sample would be a pill or tablet of known composition. 
When a tablet of equal composition was position in the other beam, the output signal would null. 
If a tablet of different composition were placed in the probe beam, a large difference signal 
would be observed. Other examples of suitable samples for comparison might include different 
regions of the ground in optical remote sensing or a piece of tissue known to be healthy in the 
case of medical diagnostics. 

In the past, the main problem with optical subtraction interferometry has been the extreme 
sensitivity to alignment of the interferometer components. Several disclosures have provided 
interferometers which are immune to small misalignments of their components, but it is thought 
that no previous disclosure has addressed the issue of providing an optical subtraction or optical 
nulling interferometer which is immune to minor misalignment of the components. The present 



disclosure extends these ideas to the problem of optical subtraction by showing arrangements of 
interferometric components which are intrinsically and optically tilt-compensated in much the 
same way as the designs disclosed by Woodruff, Turner, and others. 

Electronic subtraction of the signals from two different detectors is analogous to direct optical 
subtraction and is a practical and useful alternative, also known for many years. Generally, the 
two outputs must be out of phase for either type of subtraction to be useful. 

Figure 1 1 shows a type of four port interferometer which is constructed from a monolithic 
beamsplitter and monolithic beam combiner elements to achieve permanent alignment. It will be 
appreciated by those skilled in the art that the path of the two beams from SI and S2 will exactly 
coincide at Dl and D2 if the two blocks are made from the same thickness of material - 
independent of the exact tilt angles of the splitter and combiner. The two input beams from SI 
and S2 pass through to the two outputs which are sampled by Dl and D2. The signals from these 
detector elements are amplified and subtracted electronically. If a subtraction succeeds in 
reducing the dynamic range of the interference signal by a factor of 10, it is roughly equivalent to 
gaining 3 bits of ADC performance. If the dynamic range can be reduced by a factor of 100, then 
slight more than 6 bits of ADC performance are gained. Currently, at audio frequencies, the best 
ADCs have about 24 bits of performance, with a dynamic range on the order of 1 million: 1 To 
gain one or two orders of magnitude in performance is very significant improvement, and will 
continue to be until the SNR performance of ADCs has reached the point that the entire useful 
dynamic range of interferometric optical signals can be recorded. 

Figures 11, 12, 13 and 14 show diagrams of four-port interferometers which are intrinsically 
immune to small misalignments of the components. The outputs may be sampled by two 
detectors and subtracted, as above, or they may be recombined, and then sampled by a single 
detector. For all three of these interferometers, to maintain optical alignment, the optical path 
difference must be scanned without introducing any tih or shear into the beams. This is not a 
trivial constraint and the ramifications will be expounded below. 

An alternative to subtracting the signals from two different detectors is to record or generate a 
facsimile of the average interferogram and subtract it electronically from each subsequent 
interferogram. At any point in time, the performance of available digital-to-analog converters is 



generally better than the performance of available analog-to-digital converters by several bits. 
Thus, the ability to playback a facsimile of an interferogram is sufficiently better than the ability 
to digitize the interferogram that subtraction of the facsimile would represent a significant 
improvement over the present state-of-the-art. Preferably, the average interferogram would be 
very accurately measured by recording a long series of interferograms, thereby reducing by 
averaging the imperfections of the analog-to-digital converter. The playback of the facsimile 
would be triggered by the fringe count of the reference laser so that it was very accurately placed 
in time and retardation relative to the signal being nulled. 

The interferometers in Figures II, 12, 13, and 14 depend on parallel reflector assemblies, or a 
related design may be constructed from the type of assembly described relative to Figure 1. 
Because no tilt or shear can be introduced into either arm of such an interferometer, a preferred 
method of refractive path difference scanning uses the change of refractive index of a gas in cell. 
A gas cell may be placed in each arm, and then the pressure in one cell may be increased while 
the pressure in the other is decreased. It is anticipated that in the relatively near future, electro- 
optic scanning of interferometer path difference will become feasible. At the time of this writing, 
the available materials and configurations can only produce very small path differences. As this 
situatioji changes in the future, it is expected that the gas cells and other optical path difference 
scanners disclosed here may be replaced with the corresponding electro-optic components. 
Figure 15 shows a suitable gas cell arrangement with a piston for changing the pressures in the 
two cells. 

It is thought that there are several gases which are appropriate for use in optical path difference 
scanner by pressure change. Xenon is attractive because it is completely inert and has a high 
refractive index which means that the pressure need not be changed as much to effect the same 
difference in path. It is thought that pressure changes of approximately 100 atmospheres with a 
path length through the gas of an inch would be useful. If the path through the gas is longer, then 
the pressure change can be reduced. Nitrogen is also an attractive gas for opd scanning because it 
is inexpensive and inert. The two gas cells in Figure 15 are connected such that amount of gas in 
either cell may be varied by actuating the piston. The two gas cells have suitable windows so that 
the beams may pass through unobstructed. By varying the pressure through the actuator 
mechanism, the optical path of an interferometer may be scanned. This approach is known for 



small path differences. The utility of it here is that the opd must be scanned without introducing 
any tilt or shear in the beams. 

Figure 16 shows an arrangement with a roof reflector supported by an air bearing and driven by a 
solenoid coil. The precision of the air bearing is required in this case because the optical 
alignment is not completely compensated as it is in the other components, or when scanning with 
the gas cell arrangement. For clarity. Figure 16 does not show two roof reflectors. The second 
would be placed behind the first and is generally indicated by the dotted line. One is used for 
each of the beams for which optical path is to be scanned. This optical arrangement called 
"optical path difference scanner" can be inserted at several points in the beam path. There are 
several types of optical path difference scanners that can be used with the basic interferometers 
of Figures 11, 12, 13, and 14. Unfortunately, the refractive scanners of the nutating disk type 
cannot be used because they cause a shift in the beam. In the interferometer of Figure 1 , the shift 
is undone by a second pass through the disk in the reverse direction. In a four-port 
interferometer, the beam cannot go back on itself because it would lead to degeneracy in the 
ports - that is, it would be a two-port interferometer. 

Other suitable path difference scanners can be constructed using combinations of mechanically 
stable components such as air bearing drives, cube comer reflectors and roof reflectors. 
Numerous other methods for scanning the opd without introducing tilt or shear may be combined 
with the devices disclosed herein. 

The principles, embodiments and modes of operation of the present inventions have been set 
forth in the foregoing provisional specification. The embodiments disclosed herein should be 
interpreted as illustrating the present invention and not as restricting it. The foregoing disclosure 
is not intended to limit the range available to a person of ordinary skill in the art in any way, but 
rather to expand the range in ways not previously considered. Numerous variations and changes 
can be made to the foregoing illustrative embodiments without departing from the scope and 
spirit of the present inventions. 



